Many-body protected entanglement generation in interacting spin systems 
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We discuss a method to achieve decoherence resistent entanglement generation in two level spin systems 
governed by gapped and multi-degenerate Hamiltonians. In such systems, while the large number of degrees 
of freedom in the ground state levels allows to create various quantum superpositions, the energy gap prevents 
decoherence. We apply the protected evolution to achieve decoherence resistent generation of many particle 
GHZ states and show it can significantly increase the sensitivity in frequency spectroscopy. We discuss how to 
engineer the desired many-body protected manifold in two specific physical systems, trapped ions and neutral 
atoms in optical lattices, and present simple expressions for the fidelity of GHZ generation under non-ideal 
conditions. 



I. INTRODUCTION 

It is well-known that entangled atomic states (e.g. so-called 
spin squeezed states) potentially allow to significantly im- 
prove resolution in Ramsey spectroscopy 111, 1211 ■ Entangled 
states are also a fundamental resource in quantum information 
and quantum computation science 13tl4j]. However, in practice 
entangled states are difficult to prepare and maintain as noise 
and decoherence rapidly collapses them into classical statis- 
tical mixtures. Thus one of the most important challenges in 
modern quantum physics is the design of robust and most- 
importantly decoherence resistant methods for entanglement 
generation. 

We have recently proposed a method jsjl that allows for 
noise resistent generation of entangled states. The method 
uses the energy gap of properly designed gapped-multi- 
degenerate Hamiltonians. While the large number of degrees 
of freedom in the ground state manifold of such systems al- 
lows to create various quantum superpositions and to exploit 
rich dynamical evolution (suitable for example for precision 
spectroscopy), the energy gap prevents decoherence as local 
excitations become energetically suppressed. A simple exam- 
ple of a many-body spin Hamiltonian that illustrates the idea 
of our scheme is a multi-spin system with isotropic ferromag- 
netic interactions. These interactions will naturally align the 
spins. While all of the spins can be rotated together around 
an arbitrary axis without cost of energy, local spin flips are 
energetically forbidden. 

This paper presents a detailed analysis of this method when 
applied to trapped ions and cold atoms in optical lattices. We 
demonstrate its applicability for decoherence resistant genera- 
tion of A^-particIe Greenberger-Horne-Zeilinger (GHZ) states 
10] and the potential of the latter to be used for Heisenberg- 
limited spectroscopy. The paper is organized as follows: In 
Sec II we review one of the standard procedures used to gen- 
erate multi-particle GHZ entangled states in ion traps and 
demonstrate the detrimental effect of phase decoherence in 
such entanglement generation schemes. In Sec III we explain 
the idea of a many-body protected manifold (MPM) and dis- 
cuss how and under what conditions it can significatively re- 
duce the effect of decoherence. In Sec. IV we show the appli- 
cability of the gap protected evolution to precision measure- 



ments and discuss the significant improvement in phase sen- 
sitivity that it might provide. In section V we elaborate on the 
physical resources required for the implementation of the gap 
protected Hamiltonian in trapped ions and discuss the advan- 
tages and disadvantages of its implementation with respect to 
standard unprotected Hamiltonians. In Sec. VI we study how 
to engineer the long range interactions required for the gap 
protected evolution in optical lattice systems interacting via 
short range interactions and discuss the effectiveness of the 
MPM for GHZ generation. In Sec. VII we analyze in such 
systems the effect of non-ideal conditions such as the mag- 
netic trapping confinement and finally conclude in Sec. VIII. 



II. MULTI-PARTICLE ENTANGLEMENT GENERATION 
A. Ideal Case 

In this section we start by reviewing a method to generate 
multi-particle entangled states in a system of N spin 1/2 atoms 
by time evolution under the so called squeezing Hamiltonian. 
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As shown in Ref. |0, [Ml the Hz Hamiltonian can be imple- 
mented in trapped ions by using the collective vibrational mo- 
tion of the ions in a linear trap driven by illuminating them 
with a laser field dim. InEq.dB we used Ji"' to denote 
the collective spin operators of the N atoms: Jq°'' = \ af, 
where a — x,y,z and af is a Pauli operator acting on the i*'' 
atom and we have identified the two relevant internal states of 
the atoms with the effective spin index a = | , j . In this 
manuscript we will use units such that = 1 and assume 
to be even. 

An appropriate basis to describe the dynamics of the sys- 
tem is the one spanned by collective pseudo-spin states de- 
noted as \,J,M,(3)z ifioll . These states satisfy the eigen- 
value relations J^°^^\J, M, l3) ^ = J{J + l)\J,M,(3)z and 

ji°''\J,M,(3)z = M\J,M,f3)z, with J = N/2,...,0 and 
— J < M < J . /? is an additional quantum number associ- 
ated with the permutation group which is required to form a 
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complete set of labels for all the 2^ possible states. 

The entanglement generation process starts by preparing 
the system at i = in a fully polarized state along the x di- 
rection, \N/2, N/2)x- As any state with J = N/2 is uniquely 
characterized by M, for denoting them we omit the additional 
f3 label. Fully polarized states along x can be written as a su- 
perposition of states with different M values along the z direc- 
tion; J^Ai C^m\N/2, M)z- During the evolution the Hamilto- 
nian imprints an dependent phase to the different compo- 
nents. As the system evolves, at first the winding of the phases 

leads to a collapse of ( ji""*). However, at time xtrev = tt all 
components rephase with opposite polarization, and a perfect 



revival of the initial state is observed with (Jx 



-N/2 



(see FiglU. Specifically, the time evolution of (Ji ) for sys 
tems with TV 3> 1 can be shown to be given by; 
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Right at time to = Uev/^ the system becomes a macroscopic 
superposition of fully polarized states along the ±x direction, 
i.e. a A^-particle GHZ state of the form 

(3) 

with and (j)± real phasesgiven by — 7r/4 and tt/A + Ntt/2. 

Recent experiments (H |211 have used this type of scheme to 
generate GHZ states in trapped ions with the aim to perform 
precision measurements of ojo, the energy splitting between 
t and I levels. Ideally the use of GHZ states should enhance 
the phase sensitivity to the fundamental Heisenberg limit ifTlll . 
However, decoherence significantly limited the applicability 
of the method. 



B. Effect of decoherence 

To understand the detrimental effect of decoherence we 
first assume that the dominant type of decoherence is single- 
particle dephasing. Such dephasing comes from processes 
that, while preserving the populations in the atomic levels, 
randomly change the phases leading to a decay of the off- 
diagonal density matrix elements. We model the phase de- 
coherence by adding to Eq. (1) the following Hamiltonian 

i 

where hi (t) are assumed to be independent stochastic Gaus- 
sian processes with zero mean and with autocorrelation func- 
tion hi{t)hj{T) = 5ijf{t — t). Here the bar denotes av- 
eraging over the different random outcomes. In what fol- 
lows we will use the property that zero mean Gaussian vari- 
ables satisfy cxp[— i j^dThij)] = exp[— r(i)], with V{t) = 



FIG. 1: (Color online) Fidelity of GHZ generation vs. with and 
without gap protection. We assume a noise with step-like spectral 
density with amplitude / and cutoff frequency ujc = X ■ The dashed 
green and solid red lines are for a protected system, A = 50, with 
/ = 0.6x and / = O.lx respectively. The dot-dashed blue and 
black dotted lines are for an unprotected system, A = 0, and same 
decoherence parameters. In the inset we show { J^"' {t))/N for = 
50. 

i/gd^i /odt2./(tl-t2)lil. 

Phase decoherence causes an exponential decay of the re- 
vival peak and and N dependent decay of the fidelity, defined 

as J^{to) = {ip§"^\p{to)\Tpx"^) (See Fig. [B- Qualitatively 
the effect of phase decoherence on the evolving state can be 
understood from the energy levels of Hz (See Fig. |2j. While 
Hz commutes with both and J^"^^, Henv only commutes 
with ji"-* does not commute with Therefore, in the 

presence of phase decoherence transitions between different 
J subspaces are allowed as long as M is conserved. As all 
the states with the same |A/| value are degenerate, there is 
no energy barrier to prevent such transitions and very quickly 
the initially populated J = N/2 manifold is depleted and the 
fidelity of generating the GHZ state significantly degraded. 
Moreover, the degradation scales exponential with increasing 
N due to the exponential scalability of the number of acces- 
sible states to which the initially J = N/2 population can be 
transferred to. 

Quantitatively the effect of decoherence can be calculated 
by using the uncoupled spin basis as it diagonalizes the total 
Hamiltonian. Each state in this basis can be label as { jn''^') ~ 
|sj, 4, . . . , s%)}, where s*^ ±1 for H and fc = 1, . . . 2^. 

If at t = the reduced density matrix of the system, p, is 
given by p = 'J2k i Pk,i (0)|n('')) (n'^'^ |, after time t 



As a consequence, {ji"\t)) = Tr[ji°^/5(t)] = 
e-r(t)(4o)(t))|p^o with (4°' (t)) |r=o the expectation 
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FIG. 2: Schematic representation of the energy levels of the x-^i"^^ Hamiltonian and the effect of the different type of noise. As states with 
different J but equal \M\ are degenerate, in the presence of phase decoherence (a^ noise) they are populated during the time evolution. ax,y 
noises couple states which differ by ±1 units of M and therefore the small energy gap between them, of the order of x, naturally protects the 
system from these type of processes. 



value in the absence of noise (Eq.(|2]i). The factor e '^^'^ 
comes from the fact that the operator ji^-* = af only 
probes one-particle coherence, i.e it only connects states with 
exactly one spin flipped. 

Assuming that at t = all atoms are polarized in the x 
direction, i.e. Pk.i{0) ~ 2^^, one can show fromEq. (|5]i that 
the fidelity is degraded to 



l,k 



l + g-r(t„) 



N 



(6) 



B. Protection against phase decohernce 
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III. PROTECTED DYNAMICS 

A. Manybody protected manifold (MPM) 

Let us now consider what happens if in addition of we 
assume that there is an isotropic infinite range feiTomagnetic 
interactions between the spins so the system Hamiltonian is 



described by the Hamiltonian He ~ 



prot 



H,, where 



prot 



(7) 



Here we have assumed that all atoms are in different orbitals 
but have enough of spatial overlap that every spin interacts 
with every other spin. 

The isotropic Hamiltonian Hprot has a ground state mani- 
fold spanned by a set of iV + 1 degenerate states. They lie on 
the surface of the Bloch sphere with maximal radius J = N/2 
and are totally symmetric, i.e. invariant with respect to par- 
ticle permutations. There is a finite energy gap Eg = \N 
that isolates the ground state manifold from the rest of the 
Hilbert space. This gap is the key for the many-body pro- 
tection against decoherence. Hereunder we will refer to the 
ground state manifold as the many-body protected manifold 
(MPM). 



FIG. 3: Schematic representation of the energy levels of the Hprot + 
xJz"'^ Hamiltonian. Hprot lifts the degeneracy of the different J 
manifolds and suppresses (in the slow noise limit) couplings between 
them. So if a t = the system is in the MPM it remains there. 



The low energy spectrum of He is shown in Fig. [3] As H^ 
commutes with Hprot, in the absence of decoherence the latter 
does not affect at all the GHZ generation dynamics, however 
in the presence of decoherence the latter does significatively 
reduces the effect of local environmental noise. The protec- 
tion can be best understood by using the basis of collective 
states. In terms of collective spin operators Henv can be writ- 
ten as: 



1 

^ k=0 



(8) 



where g'^{t) 



1 



and 



?(fe) 



i " ■ Note that allowed transitions must conserve 

M as both the system and noise Hamiltonian commute with 
Jz^\ In the presence of a large energy gap Eg, one can distin- 
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guish two different type of processes: (i) Decoherence effects 
that take place within the MPM due to the collective dynam- 
ics induced by the fc = component of Henv, and (ii) tran- 
sitions across the gap induced by the inhomogeneous terms. 
The later couple the MPM with the rest of the system, how- 
ever they are nonenergy conserving process and consequently 
perturbatively weak. 

Using a perturbative analysis, and assuming that at t = 
the system lies within the MPM, the evolution of the 
projection of the density matrix on the MPM: p^/j^/ = 
z{N/2, M\p\N/2, M)z can be written as 



r)e~3(TM-7Af)^ 

(9) 



Here 



N - N M r n 

Ouit) ^ {-,M\ drH,,Ur)\^,M) ^ ^ 9\r), 

(10) 

accounts for the dynamics induced by the noise within the 
MPM and 



T^Ar/o a "'0 



M piTU.j,i3 I 2 



(11) 



takes into account the depletion of the J = N/2 levels due 

N_ 

2 



to transition matrix elements between and states 

outside the MPM: 7W^^ = ,{f,M\Henv\J,M,P),. ujj^p 
are the respective energy splittings. Because Henv is a vec- 
tor operator, according to the Wigner-Eckart theorem, Henv 
only couples the states in the MPM with states which have 
J = N/2 — 1 and thus with excitation energy AA^. 

Assuming the power spectrum of the noise, f{i^) = 
/ dte~^'^* f{t), to have a cut-off frequency Uc (e.g. f{uj) — f 
for uj < ujc and otherwise), we find that 



N^ - UP 
N 



/ 



duj 



sin(f(t< 



-AiV)/2) y 



(12) 

In the limit when the noise is sufficiently slow , i.e. 
iiJc ^ Eg, then ^^^{t) is bounded for all times, 7^^(i) < 
Jfi^'^ ^ '■^^ atomic population within the 

ground state manifold is fully preserved i.e. ^^\t) w in Eq. 

Consequently, in the slow noise limit type (ii) processes are 
energetically forbidden and only type (i) processes are effec- 
tive and therefore the noise acts just as a uniform random mag- 

netic field: if at i = p = Ea/.a/ Pm.a/WIt - ^^>(f ' ^^1' 
then after time t each component p^^^ acquires an additional 

random phase e'^^^^^*)"^*^*)) and on average 



tXt{M^-M)\-rit) 



(13) 



The factor of \'N in the denominator of 6j^f is fundamen- 
tal for the reduction of the effect of decoherence within 



the MPM. For example, it makes ji^J to decay A'' times 

slower than in the unprotected system: i.e. {Jx^y{t)) = 
e-r(0/^(j(0)(t))|,^„. 

Assuming all atoms are initially polarized in the x direction, 

PmjA^) = ^'^\/{m+n/2){m+n/2)' using Eqs. (Ql, the 

approximation (^/_|'^/2) ~ (t^)^^^^"^' '^alid in the large 
N limit and replacing the sums over M and M by integrals 
the fidelity at a given time t can be shown to be given by: 



1 



(14) 



The insensitivity of J-{t) on A^, and the A^ times slower decay 

rate of {Ji^"^) demonstrate the usefulness of MPM to generate 
a large number of entangled particles. 



C. Protection against arbitrary noise 

We now discuss the protection against spin flips, which 
can be modeled by terms proportional to erf, af in the noise 
Hamiltonian, Eq.®. First of all note that as the f and j states 
have a finite energy splitting ujq, low frequency noise associ- 
ated with such terms will be suppressed due it. However, most 
of the spin flips are generally induced by imperfections in the 
laser fields and therefore are at frequencies close to ojq, i.e. 
they correspond to low frequency noise in the rotating frame 
of the laser In the case involving GHZ state generation, the 
finite energy cost imposed by xJ? between levels with dif- 
ferent value tends to inhibit these processes as illustrated 
in Fig. 2. If in addition Hprot is present this natural pro- 
tection can be enhanced due to the fact that the energy gap 
suppresses the component of the noise that cause transitions 
between the MPM and other manifolds. More precisely, noise 
modeled as J^i h^&f, when projected into the MPM reduces 
to ^.g"(0)J° with a = x,y,z. 

In Fig. |4]we quantify the protection provided by the MPM. 
In the absence of any protection we find the finite energy cost 
imposed by H^ helps to protect the system against transver- 
sal noise. Instead of the exponential decay ( ^ q-N'^M ) 
of the fidelity observed when dephasing is present spin flips 
degrades the fidelity as ~ e-'^i^o)N/4 ^ \Q^st for the moder- 
ated A^ < 12 we have to restrict our simulations). With pro- 
tection the fidelity scales even better as the transversal noise 
is restricted to act only within the MPM . Instead of the ex- 
ponential decay of the GHZ generation fidelity with A^, with 
protection it decays as ~ p-Ar{to)N ^^^^ ^ numerical 
constant, A ^ 1/3 (for this result we do not have to restrict 
to moderated N). Hence, with Hprot we gain a factor of order 
~ V^TV. The reason why and noise degrade stronger 
the fidelity than noise (which leads just to a A^ indepen- 
dent fidelity) is that the former do not commute with H^ and 
mix states with different M quantum number Additionally in 
the figure we show that for noise with long correlation time 
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FIG. 4: (Color Online) Fidelity to create GHZ state as a function 
of N for systems in the presence of Oy noise : without MPM (green 
dashed line), without MPM but with spin echo (dot-dashed blue line), 
with MPM (red solid Une) and with both MPM and echo (dotted 
black line). We also show the degradation caused by pure dephas- 
ing in an unprotected system (long dashed purple line) for compar- 
ison purposes. For simplicity we assumed infinite correlation time: 
,f{u)) = T(5((i)/x),T = O.lx- The echo technique consisted of a 
perfect sudden vr pulse around the x direction at xt = /'^- Because 
the y components of the noise do not commute with jf the dy- 
namics was solved numerically. For the unprotected system all the 
2^ states had to be considered and we had to limit the particle num- 
ber hy N = 10. On the other hand the restriction of the dynamics to 
the MPM in the protected scheme allowed us to extend the calcula- 
tion to larger A'^ values. For the ay noise the fidelity with protection 
does not become A'^ independent, but instead scales as e^""*''^ 
(see fitted red dots). Nevertheless we do gain a factor of order ^/N 
with respect to the unprotected system. The plot also shows that the 
combination of MPM with spin echo provides the best protection. 



the use of spin echo techniques can help to further reduce the 
effect of decoherence. 



APPLICATIONS TO PRECISION MEASUREMENTS 
USING TRAPPED IONS 



IV. 



Recent experiments lll|,|2|] have generated GHZ states made 
of up to six beryllium ions and used them to perform preci- 
sion measurements of luq. For the ideal GHZ state prepara- 
tion, the spectroscopy should lead to Heisenberg-limited res- 
olution, \Sujo\ oc N^^ ifTTIl . However, in practice, even for six 
ions, the phase accuracy was significantly degraded by deco- 
herence. 

The spectroscopy lU, 01 was realized by first creating 
the desire GHZ state by applying to the initial polarized 
state, I J = N/2, N/2) z the unitary gate operation Un ^ 
gi7r/2j(o)gi7r/2j('"^g-i7r/2j(o)^ j^^^ qjj2 State was let to 
to freely precess in the z direction for time t so each atom 
accumulated a phase difference <f> = {uj — ujQ)t (in a ref- 
erence frame rotating with the frequency lu, the frequency 
of the applied field). The phase difference was then de- 
coded by measuring the collapse probability into the states 



I J = N/2, N/2)z or | J = N/2, -iV/2)^ after applying the 
unitary transformation erf. 

This generalized Ramsey sequence can be quantitatively 
described as a measure of the expectation value of the fol- 
lowing operator, O: 

(O) = (^-(^0)1 X{[cos{4>)ot - sin(0)af]|V'(io)) (15) 



with |^(to)> = e-*^^i \N/2,N/2),. 

The phase sensitivity (A^O) achievable by repeating 
the above scheme during total time T is related to the 
signal variance (A^O) = (O^) — (O)^ and given by: 



|(5wo| 



tT 



{S{0)/6<pr 



1211 . Because = 1, we just 



have to calculate (O) ~ Ti[p{to)0] to evaluate Sujq. 

Experimentally, magnetic field noise is one of the sources 
of phase decoherence. Assuming that such dephasing mainly 
takes place during the GHZ generation, as during the Ramsey 
interrogation time the atoms are essentially freely evolving, 
using Eq.(|5]l one can show that for the unprotected system 

id) = (^P^^l U^cosi^^^t - e-r(*o) sin(</.)af]|V.f 

i 

^i[^AL^ '\ >-r + h.c (16) 

Consequently the maximal phase resolution, achieved at 
4'opt = (for integer n), can be shown to be given by 



\ScJo\opt = \ScJo\sh/G, 



(17) 



with G = y/HN - l)e- 



-2r(tn 



1) and \Sujo\ 



the shot noise resolution. The factor G explains the strong 
limitations introduced by decoherence. If the time required to 
generate the GHZ state is such that G ^ I (i.e when r(to) > 
ln[\/]V]), the phase accuracy is reduced to the classical shot 
noise resolution. 

However, if instead Hprot + H, is used for the GHZ gen- 
eration, G is replaced by ^{{N - l)e-2r(to)/w + i). Due 
to the N times slower decay rate of the atomic coherences, 
the same preparation time that leads to shot noise resolution 
without protection, can lead to a sensitivity at the level of the 
fundamental Heisenberg limit with protection. 



V. IMPLEMENTATION OF THE GAP PROTECTED 
HAMILTONIAN IN TRAPPED IONS 

We now proceed to review and complement the implemen- 
tation of the protected Hamiltonian, x{J^^^'^ ^ j'z'^'^), pro- 
posed in Ref . uM . Consider a linear trap with a string of ions 
with two relevant internal levels. The ions are assumed to be 
cooled such that only the in-phase collective center of mass 
oscillation of all ions is excited. The corresponding oscilla- 
tion frequency is denoted by v. The two internal levels of 
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the ions are coupled by a laser field with a slowly varying 
Rabi frequencies and with frequency uji = ujo — being S 
the detunning from resonance. Assuming that the field couple 
all ions in the same way, we can describe the system by the 
Hamiltonian H = Hq + Hin, where ~ vo}a + 
a being the annihilation operator of the quantized oscillation 
mode. The interaction Hamiltonian is given by 

iiiu = f7j|"'e*'^*e*''('^^"""+^'=""") + h.c. (18) 

where r\ is the Lamb-Dicke parameter. The detuning 8 is as- 
sumed to be large compared to the linewidth of the resonance 
but sufficiently different from the frequency of the center os 
mass oscillation. As a consequence, the dominant processes 
are two-photon transitions leading to a simultaneous excita- 
tion of pairs of ions. 

We first assume that the ion trap is in the Lamb-Dicke limit, 
i.e., that the ions are cooled sufficiently enough, such that 
for all relevant excitation numbers n of the trap oscillation 
(n + 1)77^ <C 1 holds. In this limit one can expand the expo- 
nent in Eq. (fTsT l to first order in 77. Confining the interest to 
time averaged dynamics over a period much longer than any 
of the oscillations present in iJi„, then the oscillatory terms 
may be neglected and we are left with a more simple effective 
Hamiltonian Ref.lil] : 



n.Si = X(i<"'' - + -^JT + A(2n+ 1) Jf) (19) 

where x — ^Jz^' ^ ~ if '^^e number of phonons 

in the vibrational mode. The first term in H^f f is the desire 
protected Hamiltonian. The second term acts as an effective 
magnetic field which can be canceled by adding an external 
magnetic field or by echo techniques. The third term comes 
the from the ac Stark shift of the atomic levels due to the laser 
fields. In contrast to the standard scheme used to create Ji"''^, 
where the n dependence exactly cancels, here it does not and 
if not corrected can certainly degrade the fidelity. The degra- 
dation can be shown to be given by: 

J'ito) = ^ P„ cxp[ ^— '-^] (20) 

n 

where P„ is the initial population of the state with n phonons. 

In order to prevent this effect one has to cool the ions to 
the ground state, P„ = (5„o, which might be feasible with the 
state of the art technology or alternatively one can use spin 
echo techniques. For example if at time to/ 2 the sign of the 
laser detuning 6 is changed, then the different components will 
rotate in the opposite direction and at to the net effect due to 
the extra second and third terms in Heff will be canceled out. 

So far we have used the Lamb-Dicke and the rotating-wave 
approximation. Now we perform a more detailed analysis of 
the validity of these approximations and estimate the effect of 
deviations from the ideal situations in an actual experiment. 
To do that we follow Ref. and change to the interaction 



picture of Heff, assuming that the undesired second and third 
terms can be canceled by the techniques described above, and 
treat the small non-ideal deviations by perturbation theory. 

• Direct coupling 

Going from Eq.lfTSll to Eq.(fT9]l the off-resonant term Hd = 
. was neglected. This term correspond to direct 
single atom spin flips without any vibrational excitation. 

Changing to the interaction picture of H^f / and using the 
fact that Hd oscillates a much higher frequency that U (t) = 
^iH^fft jjj^j jjjg latter can be treated as constant in the inte- 
grals used in the Dyson series, one can show that 

T{to) = 1 - ^ sin^iSto) + msm^{6to/2) + . . . ) 

(21) 

The degradation of fidelity is a factor of N larger than the 
degradation caused by direct couplings in the standard real- 
ization of Ji"^^ where J^{to) = 1 - ^ sm^{Sto). There- 
fore it is important for the implementation of the protected 
hamiltonian to use weak laser power or to control the system 
parameter such that Sto = 2Ktt with K an integer. 

• Lamb-Dicke approximation 

In Ref.lUt] it has been shown that relaxing the Lamb-Dicke 
approximation and including higher order terms results in 
an effective Xn which depends on the vibrational number of 
phonos in the collective mode: Xn = x[l ^ if{2n + 1) + 
?7^(5/4n2 ^l^iri + 1/2)]. As this effect is global, the gap 
does not protect against it and it leads to a degradation of the 
fidelity given by 

n 

• Other vibrational modes 

With N ions in the trap, assuming that the transversal potential 
is strong enough to frozen the transversal degree of freedom, 
only the N longitudinal vibrational modes are relevant. So far 
we have assumed that only the collective center of mass mo- 
tion is excited and neglected other modes. If we include the 
effect of other modes into account, the fidelity is decreased. 
The main sources of decoherence are :a) off resonant direct 
couplings to other modes and b) reduction of the coupling to 
the center of mass mode, x, due to the vibration of the other 
modes. However all these effects are local and the gap ener- 
getically suppresses them. 

• Spontaneous emission 

Additionally a more fundamental source of decoherence 
arises from spontaneous emission effects. In typical ion trap 
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hermite conjugates and then use the projected equations of 



FIG. 5: Raman transition to a third level with atomic decay. Here A 
is the detuning of the laser fields with frequencies ujei,2 from the one 
photon resonance ,A = Uei — vi = uje2 — vi , and fii,2 are laser 
Rabi frequencies. 



experiments the | f) and | j) levels are coupled through Ra- 
man transitions to a third excited level |e) (see Fig|5]l. As- 
suming two photon resonance conditions, that is ojei — 0Je2 = 
v I — 1^2 ~ 'j-'o and A — LOei — vi = u>e2 — ^2 where A is 
the detuning of the fields from the one photon resonance and 
tJei.2 and r2i,2 are laser frequencies and Rabi frequencies re- 
spectively, the Hamiltonian of the system in the appropriate 
rotating frame can be written as 

iis = -Aj2(^le+His (23) 

i 

Hi, = n,Y.(^l+al^)+n2j2i^e + ^h) (24) 

i i 

where ct^^ = |e)ii(e|, (T|g = | i)ii{e\ and a]-^ = \ t>ii(e|. 

The decoherence processes due to spontaneous emis- 
sion can be described by means of Heisenberg-Langevin 
equations iITtIi given by: 



7-' - - 



(25) 



= (i^ii 

-(iA + re/2)a|, 
+ (ir!i - fpaj^ (26) 

(iA - r,/2)al^ + {in, + ffj(aL - a^) (27) 



-{in2 + i\l)a\^ 



(28) 



where Fg = 71 + 72 with 71 and 72 are decay rates form 
|e) to I I) and | |) respectively and the noise operators / 

have zero mean and are S correlated UtIi : ( f-i^ (t) f^^ (t')) = 



7i<5(t - and {fX{t)ffjt')) 



j2S{t ~ t')5,,k. 



In the large photon detuning limit A ^ ^1^2, 7i, 72, one 
can adiabatically eliminate the operators fr^j^ and (T^i and their 



motion to solve for (t|^. 



- 



FH + i[H, 

2 ^ 

3 



(29) 

hUt)apO) 



where FH accounts for the the Hamiltonian part of the dy- 
namics and with = ^{fH - ffj - ^{ff^ - /^J, 

= -&(/fI - He) - '^iffe - He) and = 

From the previous expressions one can estimate the degra- 
dation of the fidelity due to dephasing ( similar degradation 
of the fidelity is caused by a; or y type of noise). In the 
adiabatic limit, i.e. A ^ riiri2, 71, 72, are indepen- 
dent stochastic white noise processes with zero mean and 
autocorrelation function hf^{t)hj^{T) — -fgSijdlt — t) with 

7sp = ^jjjJ consequently the gap can not protect 

against this broad band noise. From Eq.© they will case a 
degradation of the GHZ fidelity of 



(31) 



In order to reduce the strong degradations due to this type 
of high-frequency decoherence processes one can increase 
the Raman-detuning [16] at the expense of slower evolution 
which in turn will make the system more susceptible to other 
kind of local noise (eg. magnetic field inhomogeneities). On 
the other hand the latter can be suppressed by the MPM. 

From this analysis we conclude that overhead of imple- 
menting J^°)^ — ji^-*^ instead of ji^-*^ is mainly the addi- 
tional echo technique required to remove the n dependence 
of Hp ff- Besides that, on average the same type of non-ideal 
disturbances are found in both Hamiltonians with the advan- 
tage of J*^*^)^ — Ji''''^ that the gap protects the system against 
those of them which are local in character. 



VI. IMPLEMENTATION IN OPTICAL LATTICES 

A. Engineering long-range interactions 

Up to now we have explored only the generation of an 
MPM via isotropic long-range interactions. In practice, how- 
ever, it is desirable to have a similar kind of protection gen- 
erated by systems with short range interactions such as those 
provided by cold atoms in optical lattices. These systems of- 
fer the possibility to dynamically change the Hamiltonian pa- 
rameters at a level unavailable in more traditional condensed 
matter systems. We now show how an MPM can be created in 
lattice systems and can be used to robustly generate N-particle 
GHZ states. 

We consider ultracold bosonic atoms with two relevant in- 
ternal states confined in a an optical lattice. We will assume 
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that the lattice is loaded with one atom per site, and again 
identify the two possible states of each site, with the effective 
spin index a =t, i respectively. For deep periodic potential 
and low temperatures, the atoms are confined to the lowest 
Bloch band and the low energy Hamiltonian is given by ifisll 

HbH = XI ^^T.^fl'Tj + \Y1 UaafLajiflaj - 1) 

{i,j)a- jo- 

+C/||no.jncrj (32) 

Here a^.j are bosonic annihilation operators of a particle 
at site j and state a, ha-.j = a^cLcrj, and the sum is 
over nearest neighbors. In Eq.(l32li the parameter t is the tun- 
neling energy between adjacent sites (which we assume spin 
independent i.e. spin independent lattices) and and Ua-.a' are 
the different on-site interaction energies which depend on the 
scattering length between the different species. Both Uaa' and 
T are functions of the lattice depth. We are interested in a 
unit filled lattice in the regime t ^ Ucra' where the system 
is deep in the Mott insulating phase lfl9[ I20I1 . In this limit, 
to zero order in t the ground state is multi-degenerate and 
corresponds to all possible spin configuration with one atom 
per site. A finite r breaks the spin degeneracy. By includ- 
ing virtual particle-hole excitations one can derive an effec- 
tive Hamiltonian that describe the spin dynamics within the 
one atom per site subspace ll2lll : 



Hiat — Hh + Hj — —A 



<i,j>,a 



<i,i> 



Here the coefficients are A = jU^i and x = t"^^- 



TT 



ii 



2[/^j^^). For simplicity we will now restrict the analy- 



sis to one dimensional systems and assume periodic boundary 
conditions. 

II H is spherically symmetric and in terms of collective spin 
operators it can be written as 



E 



H 



1^ f(0)2 _ 1^ ^ 

' N"^ N ^ 

k=l...N~ 



(34) 

All the iV+ 1 fully symmetric states with J = N/2 are degen- 
erate and span the ground state of Hh- Hi is not spherically 
symmetric but we can also write it in terms of collective oper- 
ators as 
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FIG. 6: (color online) Fidelity to generate a GHZ state vs A/x- In the 
inset we show {Jx'^\t)). The blue dot-dashed, dotted black, dashed 
green and solid red correspond to A = 0, 5, 10, 20 respectively. The 
plots are obtained by numerical evolution of Eq. l l33t for A'^ — 10. 



of Hi on it, which corresponds to VHi = Xejf^'^ — 
with Xe = i^ZT' is effective and Hi acts as a long range 
Hamiltonian. Here we used the relation 7'fe^o['^i'^'''^2 '^''l = 
+ 4(^^, with V the projection into the J = N/2 
subspace. The non zero projection of the latter term comes 
from the fact that the operators Ji'^^ Ji and jf^"^ are not 

independent as they satisfy the constrain J^fc^Jo^ j'^^^ '^^ = 
iVV4. 

In Fig. |6]we contrast the dynamical evolution of a system 
in the presence and absence of Hh assuming at time t ^ 
all the spins are polarized in the x direction. If only the Ising 
term is present, A 0, it induces local phase fluctuations that 
leads to fast oscillations in (Ji^^) = N/2cos^[2xt]. On the 
other hand, as the ratio A/x increases, the isotropic interac- 
tion inhibits the fast oscillatory dynamics and instead (ji^"*) 
exhibits slow collapses and revivals. For A ^ x the dynamics 
exactly resembles the one induced by Hz and at Xet = tt/2 
the initial coherent state is squeezed into a GHZ state. 



B. MPM in lattice systems 



Hr 



4x 

N' 



4x 
N 



k=l...N-l 



2'Kk 



(35) 

If the condition x ^ is satisfied, which can be engineered 
in this atomic systems by means of a Feshbach resonance, the 
effect of the Ising term can be studied by means of perturba- 
tion theory. Assuming that at t = the initial state is pre- 
pared within the J — N/2 manifold, a perturbative analysis 
predicts that for times t such that xt < A/x, Hh confines the 
dynamics to the ground state manifold and transitions outside 
it can be neglected. As a consequence, only the projection 



Hh also provides protection against phase decoherence. 
However, Hh is not as effective as Hprot because the en- 
ergy gap between the MPM and the excited states of Hh van- 
ishes in the thermodynamic limit as Eg X/N^ . This is 
a drawback of the short range Hamiltonian for the purpose 
of fully protecting the ground states from long wave length 
excitations. Note however that one dimensional systems are 
the worst scenario as for higher dimensions the gap vanishes 
as iV^^/'' with d the dimensionality of the system. Never- 
theless, the many body interactions can still eliminate short- 
wavelength excitations since in the large N limit they remain 
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separated by a finite energy gap, 8A. 



We quantify the effectiveness of the MPM to protect the 
system against Henv by using time dependent perturbation 
theory. For this analysis we restrict to the limit A 3> x 
where the Ising term can be treated as an effective XeJz^^^ — 
Hamiltonian. In this limit a convenient basis to study 
the quantum dynamics is the collective spin basis. Assum- 
ing that at f = the system lies within the J = N/2 
manifold, the evolution of the matrix elements p^^jq^ = 
;,{N/2, M\p\N/2, M)^ can be written as 



(38) 



In Figl?] we plot T calculated from Eq.dJSlias a function of 
N. In the lattice the fidelity is degraded as N grows because 
the gap decreases and the generation time increases with A^. 
Moreover, an abrupt drop of the fidelity occurs at the value of 
N at which ~ lo,.. 



(36) 

where the random phase, given by Eq.JTOb. characterizes the 
dynamics induced by the noise within the MPM and 7^^^ (t) = 

f3 I lo dTM^/pe'-'^'^^"-''f> p, takes into account the de- 
pletion of the J ~ N/2 levels due to transition matrix ele- 
ments with states outside the symmetric manifold: M.^/^ = 

z{^tM\ Henv I J, M, fi) z. w'"^ are the respective energy split- 
tings. Up to this point the expressions are structurally identi- 
cal to the ones obtained for long range interactions. The dif- 
ference appears in the evaluation of ^lat- In contrast to Hprot, 
not only the excitation frequencies '"g are not degenerated 
but also they become smaller as N is increased. As a conse- 
quence Eq. (fT2] | is replaced by the following equation for the 
lattice system 



N^ - 4M 



N{N -1) 



fe=i 



dcu 



sm{t{uj - AE'')/2) 



VII. NOISE AND DECOHERENCE IN LATTICE SYSTEMS 

In the previous section we used the effective Hamiltonian 
given by Eq. ( |33] | to study the GHZ generation in lattice sys- 
tems. Here we perform a more detailed analysis of its validity 
and estimate the effect of deviations from the ideal situations 
in an actual experiment. For this analysis we restrict to the 
limit A ^ X where the Ising term can be treated as an effec- 
tive XeJz^^^ — Hamiltonian. 



with E'^ the excitation energies of the states that belong to the 
J = N/2 - 1 manifold given by AE'' = 8Xsin^{T:k/N), 
with fc = 1, . . . , — 1 I22I1 . From Eq.(l37]i we can estimate 
the degradation of the fidelity due to phase decoherence as: 




A. Particle-hole excitations 

Deriving Eq. ( |33] | from the Bose-Hubbard Hamiltonian we 
2 only included virtual-particle hole excitation. However, dur- 
(37)ing the time evolution real transitions from singly to doubly 
occupied states can take place and they degrade the fidelity. 

To account for these effects, we write the manybody wave 
functionas = J^n CnlV'n>+Em 5m|0m), wherc IV-ri) 

span the Hilbert space with one atom per site and | (/),„) span 
the subspace with one particle and one hole adjacent to each 
other and N — 2 singly occupied sites. The latter are the states 
that directly coupled to \ipn) through tunneling. Solving the 
time dependent Schrodinger equation from the Bose Hubbard 
Hamiltonia, using the assumption that Ua,cr' ~ U and that 
at time i = no doubly occupied states are populated one 
obtains 



25 50 75 100 125 150 

N 



(39) 



Here we also assumed that {C„} change at a rate much 
smaller than Ua,a' and treated them as constants during the 
time integration. Eq.(|39]) yields the following lost of fidelity 
due to real particle hole excitations: 



FIG. 7: In the presence of phase decoherence the fidelity of the GHZ 
state preparation in the lattice is degraded as N grows because in the 
lattice the gap decreases and the generation time increases with A'^. 
In this plot we assumed the limit x ^ A where Eq. B8b holds and 
used a system with cjc = ^ = lOOx and T = O.Olx- At A = 90, 
AEg = and it explains the drop of for N > 90. 



J^{to)-l - ^sm\Uto/2). (40) 

remembering that Xeto ~ 7r/2. As long as x/U ^ 1, we con- 
clude that particle-hole excitations do not significantly affect 
the GHZ generation. 
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B. Magnetic confinement 

In Eq. (|32] | we assumed a translationally invariant system. 
However in most of the experiments an additional quadratic 
magnetic confinement is used to collect the atoms. Actually is 
due to this quadratic potential that a unit filled Mott insulator 
has been experimentally realized. In its absence it would be 
difficult to create a unit filled Mott insulator as in an homo- 
geneous system it only takes place when the number of atoms 
is exactly equal to the number of lattice sites. A drawback 
of the magnetic confinement is that it generates always super- 
fluid regions at the edge of the cloud, so only a fraction of 
the total trapped atoms located at the trap center has to be se- 
lected as the quantum register. Assuming we work on this unit 
filled Mott Insulator subspace, here we quantify the effect of 
the magnetic potential in the GHZ generation in the A/x ^ 1 
limit. 

The magnetic confinement is accounted for by adding a 
term W J^j a J^^-crj Bose Hubbard Hamiltonian. W = 
\/2mLd'^a\ with m the atom mass, lot the frequency of the 
external trapping potential and a/, the lattice spacing. This 
term modifies the global coupling constants A and x when 
the effective Hamiltonian is derived and make them site de- 



pendent, A — > A J = T . 1 1 and x 



-1 
•j.TT 



U, 



2u; 



..,^^^). HereC/.,.., U,,. / {U^^, ^W\2i + if ). 
Assuming that the gradient of the external potential is weak 
compared to the on site interaction energy, as is in general the 
case for current experiments, the effective Hamiltonian in the 
presence of the magnetic trap becomes 



Hi = 



T, 



Hlat + 

<ij> 

[73 



(41) 
(42) 

(43) 



The corrections on the fidelity of the GHZ state introduced 
by H]^ can be estimated by calculating the effective projection 
of it on the MPM: VHf . As the latter is just proportional 
to the identity matrix /, VH^ = J2i ^j^' i'^ effects is just a 
global phase and it does not cause any main degradation of 
the fidelity. Similarly any other perturbation induced by local 
fluctuations in the magnetic field or the lasers used to generate 
the lattice become irrelevant thanks to the MPM. 

From this analysis we conclude that except from sponta- 
neous emission or heating mechanisms the MPM effectively 
protects lattice systems against common non-ideal situations 



encountered during their experimental realization. On the 
other hand lattice-based GHZ state generation faces the scal- 
ability problem due to the fact that the gap decreases and the 
generation time increases with increasing N. 

VIII. CONCLUSIONS 

We have in this paper evaluated the possibility for a robust 
preparation of multi-particle GHZ entangled states of trapped 
ions or cold atoms in optical lattice by generating a decoher- 
ence free multi-level manifold corresponding to the ground 
levels of properly designed Hamiltonians. The MPM is iso- 
lated from the rest of the Hilbert space by an energy gap which 
energetically suppresses any local decoherence processes. We 
have presented analytical estimates for the fidelity of the GHZ 
preparation. 

In trapped ions we demonstrated that the fidelity can be 
significantly better than the one achievable without any gap 
protection and therefore that our scheme is in the position to 
improve the spectroscopy resolution in current Ramsey spec- 
troscopy experiments. 

We also showed that cold atoms in optical lattices inter- 
acting via short range interactions can be utilized to engineer 
long range interactions which in turn can be used for gener- 
ating many-body entanglement. We calculated the effects of 
non-ideal conditions and concluded that the main restriction 
in these systems is the scalability as the MPM protection de- 
grades with increasing N. 

The scalability certainly limits the use of lattice systems for 
massive entanglement generation, however it is not a prob- 
lem for recent quasi-one dimensional experiments 1230 where 
an array of ID tubes with an average of 18 atoms per tube 
has been realized. In such systems therefore it should be 
possible to create few-particle collective entangled states us- 
ing our scheme and to perform proof-of -principle experiments 
demonstrating the improvement of spectroscopic sensitivity. 

We emphasize that, even though we have limited the dis- 
cussion to ensembles of spin = 1/2 particles, the MPM 
ideas can be straightforwardly generalized to systems com- 
posed of higher spin atoms. Besides entanglement genera- 
tion, the MPM might have also important applications for the 
implementation of good storage memories using for example 
nuclear spin ensembles in solid state ll24ll or photons ll25ll . 
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